Disorder Driven Melting of the Vortex Line Lattice 
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We use Monte Carlo simulations of the 3D uniformly frustrated XY model, with uncorrelated 
quenched randomness in the in-plane couplings, to model the effect of random point pins on the 
vortex line phases of a type II superconductor. We map out the phase diagram as a function of tem- 
perature T and randomness strength p for fixed applied magnetic field. We find that, as p increases 
to a critical value p c , the first order vortex lattice melting line turns parallel to the T axis, and 
continues smoothly down to low temperature, rather than ending at a critical point. The entropy 
jump across this line at p c vanishes, but the transition remains first order. Above this disorder 
driven transition line, we find that the helicity modulus parallel to the applied field vanishes, and 
so no true phase coherent vortex glass exists. 
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p and temperature T. Increasing p at fixed B is believed 
to play a similar role as the more physical case of increas- 
ing B at fixed p. Taking great care to achieve proper 
equilibration, and measuring thermodynamic derivatives 
of the free energy, we find in our model a single first order 
phase boundary. At small p, the transition, T c (j>), is a 
thermally driven melting of the vortex lattice. Increasing 
p, there is a maximum critical p c above which disorder 
destroys the vortex lattice; when p = p c , the transition 
line T c (p) turns parallel to the T axis and smoothly con- 
tinues down to low T. No critical point is found to sepa- 
rate the thermal from the disorder driven sections of the 
phase boundary. For p > p c we find no evidence for a 
true phase coherent vortex glass. 

The model we study is given by the Hamiltonian 



Experimental |]lj-^f, theoretical |^-|6| and numerical 
studies have argued that the effect of intrinsic point 
impurities on otherwise clean single crystal samples of 
high T c superconductors leads to a H — T phase diagram 
with the following generic form. At low magnetic fields 
H, an elastically distorted vortex lattice (the "Bragg 
glass" (J]) undergoes a first order melting transition to a 
vortex liquid as temperature T is increased. This melting 
line T C (H) continues as H is increased, until an "upper 
critical point," T ucp , is reached above which sharp discon- 
tinuities in measured quantities become smeared. Upon 
increasing H at lower temperatures, T < T ucp , the vortex 
lattice undergoes a transformation to a disordered vortex 
state along a line, H sp (T), characterized by the "second 
magnetization peak" where critical currents show a sharp 
increase. As T increases, the H sp (T) line continues to the 
vicinity of T ucp . In BSCCO, H sp (T) is only weakly de- 
pendent on T Very recently, several experiments || 
have provided strong evidence that the H sp (T) line in 
BSCCO is associated with a thermodynamic first order 
phase transition. Whether the disordered state above 
H sp (T) is a "vortex glass" M, characterized by true su- 
perconducting phase coherence and separated from the 
vortex liquid by a sharp phase transition, or whether it 
is a dynamically frozen state that smoothly crosses over 
to the vortex liquid, remains a topic of controversy E|. 

Since many of the experimental and numerical stud- 
ies focus on dynamical probes, from which it can some- 
times be difficult to infer a true equilibrium phase tran- 
sition, and analytical models must resort to Lindcmann 
or other simplifying approximations, it is important to 
establish the true equilibrium phase diagram within a re- 
alistic model system. Towards this end we have carried 
out extensive Monte Carlo (MC) studies of the uniformly 
frustrated three dimensional (3D) XY model pH , with 
uncorrelated quenched random couplings to model point 
impurities. Applying a fixed magnetic field B we map 
out the phase diagram as a function of disorder strength 
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where 9i is the phase of the superconducting wavefunc- 
tion on site i of a 3D periodic cubic grid of sites, the 
sum is over all bonds in directions p — x.y, z, and 
Ai^ — (27r/0o) /* A ■ is the integral of the mag- 
netic vector potential across bond ifi, where V x A = Bz 
is a fixed uniform magnetic field in the z direction. To 
model uncorrelated random point vortex pinning in the 
xy planes, we take, 
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The coupling between planes J z is uniform, while each 
bond in the xy plane is randomly perturbed about the 
constant value J± ; the are independent Gaussian ran- 
dom variables with (e 1M ) = 0, (ef„) = 1. The parameter p 
controls the strength of the disorder. For computational 
convenience we choose J z /J± — 1/40, with a vortex line 
density of / = o^B/^q = 1 /20, where a± is the grid spac- 
ing in the xy plane, and 4>o — hc/2e is the flux quantum. 
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Our system size is L x — L y = 40, with L z = 16. To check 
finite size effects, we have also considered L z — 24 and 
32 for certain cases. Our runs are typically 1 — 10 x 10 7 
MC sweeps through the entire lattice near transitions. 
Our results below are for a single realization of the dis- 
order only. The extremely time consuming nature of our 
simulations excluded any serious attempt at disorder av- 
eraging. We have, however, carried out similar analyses 
for two other realizations of the disorder and have found 
qualitatively the same behavior. 



uid. For the stronger p = 0.09, the disorder sufficiently 
breaks the degeneracy of the two lattice orientations, so 
that we see only coexistence between one particular lat- 
tice orientation and the liquid. The repeated hopping 
between lattice and liquid in Fig. 1 verifies that we are 
well equilibrated. In this manner, by varying p, we deter- 
mine the melting phase boundary T c (p) shown in Fig. 2. 
We include in Fig. 2 several data points for systems with 
L z = 24, showing only a small shift in the phase bound- 
ary as L z increased. 



0.25 
0.20 
U7 0.15 
£ 0.10 
0.05 
0.00 



(a) 



p=0.06 



10 15 20 25 

MC sweeps t ( x 10 6 ) 



30 35 



p=0.09 j 




20 30 40 50 60 
MC sweeps t ( X 10 6 ) 

FIG. 1. Si (•) and S[ (o) vs. number of Monte Carlo 
sweeps t for disorder strengths (a) p = 0.06, T c = 0.255, and 
(b) p = 0.09, T c = 0.200. The system width is L z = 16. Each 
data point is an average over 2 16 sweeps. 

In the pure model, p — 0, the low temperature vortex 
line lattice has a first order melting transition to a vortex 
line liquid Jl2] ]. To map out this melting transition line 
in the p — T plane we fix the disorder strength p, and cool 
down from high T, until we reach a temperature T c (p) 
at which we observe repeated hopping back and forth 
between coexisting vortex liquid and lattice phases. To 
detect the vortex lattice, we measure the in-plane vortex 
structure function (normalized for convenience so that 
5(0) = 1), 
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where n z is the vorticity in the xy plane, and N = 
L x L y L z . S(kjJ will have Bragg peaks at the reciprocal 
lattice vectors {K} of the vortex lattice. We find that the 
vortex lattice always orders into the same periodicity as 
that of the pure p = case, where there are two possible 
lattice orientations related by a 90° rotation. Defining 
Si as the average of S(K) over the six smallest non-zero 
{K} for one lattice orientation, and S[ as that for the 
other orientation, we identify the vortex lattice as states 
in which either Si or S[ is large, according to which of 
the two lattice orientations has formed; in contrast, in 
the vortex liquid, both Si and S[ are small. In Fig. 1 we 
plot Si and S[ vs. MC simulation time t at T c (p) for two 
disorder strengths. For the weaker disorder, p = 0.06, we 
see both orientations of lattice coexisting with the liq- 
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FIG. 2. Vortex lattice melting phase boundary in the p — T 
plane, for L z = 16 (o) and L z — 24 (•). 

We now determine if the melting transition T c (p) re- 
mains 1st order, asp increases. First order transitions are 
characterized by discontinuous jumps in thermodynamic 
quantities. Here we consider the average energy per site 
E (the conjugate variable to temperature T), and a vari- 
able Q defined to be conjugate to the disorder strength 
P, 
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where F is the total free energy. 

To see if there is a discrete jump in E or Q at T c (p), 
we use the values of Si and S[ to sort microscopic states 
as either vortex lattice or liquid. We can then compute 
the properties of each phase separately. In Figs. 3a and 
3b we show semilog plots of the histograms P(ASi) of 
values of A Si = Si — S[ encountered during our simu- 
lation at T c (p), for the two cases of Fig. 1. In Fig. 3a, 
for p — 0.06, we see separated peaks for the liquid at 
A Si = 0, and for the two lattice orientations at finite 
positive and negative values of ASi . In Fig. 3b, for the 
stronger p = 0.09, we see only peaks for the liquid and 
one of the two lattice orientations. Fitting these peaks 
to empirical forms (Gaussian for the lattice, exponential 
for the liquid; these are the solid lines in Figs. 3a,b), we 
determine the relative probability for a state with a given 
value of ASi to belong to the liquid phase, or either of 
the two orientations of the lattice phases. Sorting though 
our microscopic states we probabilistically assign each to 
one of these three phases. We then plot the histograms 
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of E and Q values separately for each phase. 
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FIG. 3. Histograms of AS\ = S\ — S[, energy E, and dis- 
order conjugate Q for p = 0.06, T c = 0.255 and p = 0.09, 
T c = 0.200, for L z = 16, at the melting temperature T c (p). 
(o) and (+) are for the two lattice orientations, (x) are for 
the liquid. 

In Fig. 3c we show the histograms P(E) for p = 0.06. 
While the two lattice orientations have the same energy 
distribution, there is a clear difference between the liquid 
and lattice. This results in a finite energy jump AE, and 
hence a finite entropy jump AE/T, between the liquid 
and lattice. The melting transition is clearly first order. 
The histograms P{Q) for p ~ 0.06 are shown in Fig. 3e, 
where we also see a finite jump AQ between liquid and 
lattice; for Q the disorder also couples differently to the 
two lattice orientations. 

In Fig. 3d we show the histograms P(E) for the more 
strongly disordered p — 0.09, where only one lattice ori- 
entation is found. In contrast to Fig. 3c, we find that 
the energy distributions of the liquid and the lattice are 
now identical! Thus there is no energy jump, and hence 
no entropy jump, and moreover no specific heat jump, in 
going from liquid to lattice. However the histograms of 
P{Q) plotted in Fig. 3f remain clearly different for liquid 
and lattice, and so there remains a finite jump AQ at 
the melting transition. The transition at p = 0.09 there- 
fore remains first order, even though the entropy jump 
has vanished. Combining AE = with the Clausius- 
Clapeyron relation, we conclude that the melting line 
must now be perfectly parallel to the temperature axis 
p3[ , and the transition becomes disorder, rather than 
thermally, driven. As shown in Fig. 2, we have been able 
to follow the melting line from where it first turns parallel 
to the T axis, down to several lower temperatures. 

In Fig. 4 we plot AE and AQ vs. p along the melting 
line. That AE and AQ never simultaneously vanish, 
indicates that no critical point exists along the melting 
line. 

Although our simulations are for a specific fixed value 
of B, if we assume that the above result continues to 
hold for general values of B, then it must be true that the 
phase diagram in the B—T plane at fixed p must similarly 
turn parallel to the T axis at low temperatures. The 
point where the melting line first turns parallel to the T 



axis has many features in common with an "upper critical 
point" ; discontinuities as T varies across the melting line 
below this point, will cease to exist as T varies above this 
point. Yet there is no true critical end point and the first 
order melting line extends continuously down to lower 
temperatures. 
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FIG. 4. Jumps AE and AQ, vs. p, along the melting 
line, T c (p), for L z — 16. Subscripts "1" and "2" refer to the 
two possible orientations of the vortex lattice (solid and open 
symbols respectively). 

Next, we investigate the superconducting phase coher- 
ence in our model, by computing the helicity modulus 
Jll| |, T 2 , parallel to the applied magnetic field. T z — 
indicates the absence of phase coherence. In Fig. 5a we 
show T 2 vs. T, for two different disorder strengths 
p < p c = 0.09, comparing systems with L z — 16 and 
L z = 24. In both cases we see a discontinuous jump in 
T 2 at the melting transition (at T c we compute T 2 sep- 
arately for the coexisting lattice and liquid phases, using 
the decomposition of our states according to ASi). As 
L z is increased, we see that Y 2 vanishes in the vortex 
liquid; the vortex lattice melting thus marks the loss of 
superconducting phase coherence. 
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FIG. 5. Longitudinal helicity modulus T z /J z vs. T for 
disorder strengths (a) p — 0.02 and p = 0.06 < p c (open 
symbols and solid lines are for L z — 16; solid symbols and 
dashed lines are for L z — 24) and (b) p — 0.12 > p c . 

We now consider p > p c to see if a phase coherent 
vortex glass state might exist above the disorder driven 
melting of the vortex lattice. In Fig. 5b we show T z vs. 
T for p = 0.12 and system sizes L z — 16 and L z — 
32. Although it becomes extremely difficult to equilibrate 
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our L 2 = 32 system at such low T, we see a dramatic 
finite size effect strongly suggesting that T z decreases to 
zero as L z increases. We therefore find no evidence for a 
true phase coherent vortex glass. Our result agrees with 
conclusions from dynamical simulations of an interacting 
vortex line model by Reichhardt et al J8|. In contrast, 
recent work by Kawamura found a finite vortex glass 
T c in numerical studies of a model similar to Eq.(Q), but 
at a higher field density / = 1/2-k and uniformly 
distributed on [0, 2]. It is not clear if this disagreement is 
due his much stronger random pin strengths, or possible 
finite size effects in his systems of L < 16. A vortex glass 
transition has also recently been found for an interacting 
vortex line model with / = 1/2 and very strong pinning 

S 

Finally, in Fig. 6 we plot the longitudinal phase angle 
correlation length, £ z , in the liquid T > T c (p), as deter- 
mined by fitting the correlation function, 



c(*) = £< 



.J) 
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to an exponential decay for z < L z /2. We see that for 
fixed T, £ z decreases as p increases. However, since T c (p) 
decreases as p increases, the value £z(T c (p)) at melting 
increases as p increases. Even in the disordered state 
above p c we find that £ z can be as large as in the liq- 
uid just above the thermally driven melting line. This 
suggests that the disorder driven melting should not be 
thought of as a layer decoupling transition. 
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FIG. 6. Longitudinal correlation length £ z vs. T in 
the vortex liquid phase for (a) several values of disorder 
p < p c = 0.09 at L z = 24; and for (b) p = 0.12 > p c at 
L z — 32, in comparison to p — 0.07 < p c from (a). 

As we were finishing this work, we learned of sim- 
ilar work by Nonomura and Hu |]l6| , using the same 
model (|l|) with a slightly different scheme for the point 
randomness and a much weaker inter-planar coupling, 
J z /J± = 1/400. Using different methods they too find 
a disorder driven 1st order vortex lattice melting line 
nearly parallel to the T axis at low T. However they also 
claim to find a 1st order "vortex slush" to liquid tran- 
sition extending to higher disorder from the thermally 
driven melting line, as well as a vortex glass to vortex 
slush transition at lower T. We too find a peak in spe- 
cific heat for p > p c that lies at a T in the vicinity of 



the thermally driven melting line, however we have in- 
terpreted this as a smooth crossover rather than a true 
phase transition. 
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